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We present results from a series of NRQCD simulations of the charmonium system, both in the quenched 
approximation and with n/ = 2 dynamical quarks. The spectra show evidence for quenching effects of ~ 10% 
in the S- and P-hyperfine splittings. We compare this with other systematic effects. Improving the NRQCD 
evolution equation altered the 5-hyperfine by as much as 20 MeV, and we estimate 0(a s ) radiative corrections 
may be as large as 40%. 



1. The current status of Charmonium 

Past lattice NRQCD simulations of charmo- 
nium have substantially underestimated the hy- 
pcrfine spin-splittings. The introduction of 0(v 6 ) 
corrections to the NRQCD Hamiltonian caused a 
dramatic decrease in the S-hyperfine splitting jl] , 
which fell below 50% of the experimental value. 
These 0(v 6 ) simulations also demonstrated a 
large dependence on the definition of the tadpole 
correction factor. Even in the less-relativistic T 
system, the same highly-improved NRQCD ac- 
tion has not provided conclusive agreement with 
experiment J|||. 

Relativistic heavy quark actions have fared 
somewhat better, though still significantly un- 
derestimate spin splittings. A tadpole- and SW- 
improved charmonium simulation by UKQCD [|J 
gave the 5-hyperfine splitting roughly 40% below 
experiment. Unquenching should increase the hy- 
perfine splittings — in [pi, quenching effects were 
projected to be as large as 40%, though 5-15% 
seems more likely ||p|. 

We describe a series of NRQCD charmonium 
simulations with both quenched and unquenched 
lattices. The effects of different tadpole correc- 
tion schemes are examined, and we find a rough 
estimate of the effect of 0(a s ) radiative correc- 
tions. Finally, we note a sizeable systematic ef- 
fect due to an instability in the usual heavy-quark 
evolution equation. 
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2. The NRQCD and Gluon Actions 

The NRQCD quark propagator is found from 
an evolution equation, 
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where the Hamiltonian H — H$ + 5H is 
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A tilde signifies the use of improved versions of 
the lattice operators jij, and all operators are 
tadpole-improved. Much evidence suggests that 
the Landau definition of the tadpole factor is su- 
perior to the plaquette definition 0,0]. We have 
used both in our simulations. 

The gluon action used for the quenched simu- 
lations is tadpole and 'rectangle' improved 0, 
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where P Aty (n) and i? M „(ri) are the traces of 1 x 1 
plaquettes and 2x1 rectangles of link operators 
respectively. We used smeared operators for the 
1 5 , 3 Si, 1 Pi, 3 Pq, 3 -Pi and 3 P 2 states, previously 
detailed in Jl]]. 

3. Simulation Details 

We performed simulations with the NRQCD 
Hamiltonian above, truncated to 0(v 4 ) and 
0(v 6 ) , with both the Landau and plaquette tad- 
pole factor. The MILC collaboration generously 
provided us with an ensemble of 200 unquenched 
gauge field configurations of volume 16 3 x 32, 
created with the Wilson action at j3 = 5.415 
and rbf = 2 staggered dynamical quarks with 
m = 0.025 (m ps /m v ~ 0.45). 

We produced ensembles of 100 quenched con- 
figurations of volume 12 3 x 24, with both Landau 
and plaquette tadpole factors. Using (3^ — 2.1 
and (3p — 2.52 gave roughly the same lattice 
spacing as the unquenched configurations. As 
heavy mesons are small, the difference in lattice 
volume between the quenched and unquenched 
configurations should not effect our results signif- 
icantly. The lattice spacing was determined for 
each ensemble using the spin-averaged P-S split- 
ting. The bare charm mass Mq was tuned by 
requiring that the kinetic mass of the -"-So char- 
monium state, defined by E(p) — E(0) — p 2 /2A4 k , 
agrees with the experimental mass of the rj c . The 
simulations parameters are shown in Table [|. 

To increase statistics we calculated eight me- 
son propagators on each configuration (four spa- 
tial sources, two time slices). Single-exponential 
fits were used for meson masses, while the S- 
hyperfine splitting was extracted from a ratio fit. 

4. Results 

Table || contains the final results (in GeV) for 
the quenched and unquenched charmonium mass 
fits. Evidently NRQCD simulations of the char- 
monium system have a number of issues yet to be 
resolved. This is most readily seen in the hyper- 
fine splittings, collected in Figures and |. 

The 0(v 6 ) corrections lead to a disturbingly 
large decrease in the S-state hyperfine splittings, 



taking them as much as 60% further away from 
the experimental values as first noted in Jj[]. 
The plaquette-tadpole results are strikingly bad, 
where the 3 P states appear in the wrong order. 
This reversal is corrected in the Landau-tadpole 
simulations, though the hyperfine splittings are 
still badly underestimated. 

The large discrepancies in spin-dependent split- 
tings would be less worrisome if quenching were 
seen to have a considerable effect on the spec- 
trum, as suggested in ||. Sadly this does not 
seem to be the case. There is some evidence for a 
~ 10% quenching effect in the 0(v 6 ) S-hyperfine 
data, though given the apparent size of other sys- 
tematic uncertainties no great significance can be 
attached to these differences. 

5. Other Systematic Errors 

These results suggest quenching errors alone 
will not account for the discrepancies in the spec- 
trum. We have seen, as others have previously 
large differences between the Uq and Uq results. 
This is not surprising, as terms in the NRQCD 
Hamiltonian linear in E or B differ by as much 
as 30% between the different tadpole schemes in 
our simulations, an effect at least as important as 
quenching effects. 

5.1. Radiative Corrections 

We expect some effect on the spectrum from 
high-momentum modes that are cut off by the 
finite lattice spacing. These effects appear as 
0(a s ) corrections to the coefficients Cj in the 
NRQCD Hamiltonian. To date, corrections have 
been found for just three 0(v 4 ) terms: the c\ 
and C5 'kinetic' terms H and the spin-dependent 
C4 term jlo|. The corrections are ~ 10% for the 
bottom quark, but rise dramatically as the bare 
quark mass falls below one (in lattice units). 

To roughly estimate the effects of 0(a s ) cor- 
rections to all coefficients, we replaced the c, = 1 
with Cj = 1 ± a s . On our lattices, a s (Tt/a) ~ 

2 

■f- ~ 0.2. We used the calculated values for ci, 
C4 and C5, and varied the remaining coefficients 
between 0.8 and 1.2. This changed the S- and 
P-hyperfine splittings by as much as 40% — while 
this is only a crude estimate, it is clear that the 
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Table 1 



Parameters used in charmonium simulations. 



/3 


u 




a (fm) 










Quenched 








2.52 


0.874 




0.168(3) 


0.81 




2.10 




0.829 


0.181(3) 


1.15 
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Unquenched 






5.415 


0.854 




n i kK'W 




u 


5.415 




0.800 


u . 1UOI O 1 


1.15 
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Table 2 












Quenched (top) and unquenched (bottom) charmonium 


masses in GeV. 






State 




u o 






Expt 




0{v 4 ) 


0(v 6 ) 


0(v 4 ) 


0(v 6 ) 




3 s 1 


3.086(2) 


3.022(2) 


3.066(2) 


3.036(1) 


3.098 


'Pi 


3.517(17) 


3.470(17) 


3.499(18) 


3.479(17) 


3.530 


3 Po 


3.439(16) 


3.522(20) 


3.426(15) 


3.441(17) 


3.400 


3 Pi 


3.486(17) 


3.488(17) 


3.483(18) 


3.478(18) 


3.520 


3 P2 


3.576(17) 


3.449(16) 


3.528(22) 


3.475(17) 


3.570 


Si— Sq 


0.106(3) 


0.042(2) 


0.086(2) 


0.056(2) 


0.118 




3.087(2) 


3.028(2) 


3.068(2) 


3.043(1) 


3.098 


'Pi 


3.514(17) 


3.475(17) 


3.500(17) 


3.486(16) 


3.530 


3 Po 


3.456(15) 


3.518(20) 


3.437(15) 


3.445(15) 


3.400 


3 P1 


3.501(17) 


3.499(17) 


3.490(17) 


3.479(17) 


3.520 


3 P 2 


3.548(21) 


3.462(16) 


3.515(21) 


3.489(17) 


3.570 




0.107(2) 


0.049(1) 


0.087(2) 


0.062(1) 


0.118 



effects of radiative corrections are easily as im- 
portant as quenching effects for heavy-quark sys- 
tems. Accurate determinations of the remaining 
0(a s ) corrections are sorely needed. 



5.2. Improving the Evolution Equation 

Equation ([!]) contains better-than-linear ap- 
proximations to the exponential e Ht for the H 
terms but only a linear approximation for the 5H 
terms. As the high-order corrections are large for 
charmonium, this approximation is perhaps too 
severe. Lewis and Woloshyn made a similar sug- 
gestion in their NRQCD study of D mesons [jul . 

We examined this possibility for the 0(v 6 ) 
terms by using an improved form for the evolu- 
tion equation that incorporates a 'stabilisation' 



parameter for the correction terms, replacing 

(l-a^(l-^)\ (6) 

A simulation using this replacement, with s$ — 
4 and all other parameters as for the previous 
Landau-tadpole quenched simulations, showed 
that the S'-hyperfine splitting increased by as 
much as 20 MeV. Statistical uncertainties in the 
P-hyperfine splittings were large, though a simi- 
lar increase seems likely. 

6. Conclusions 

While NRQCD is considered quite unreliable 
for the charm quark, the variety of relativis- 
tic quark actions also tend to underestimate the 
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Figure 1. S'-hyperfine splitting. 4, 6 = 0(v 4 ) 
0(v 6 ) ; P, L = Plaquctte, Landau tadpoles. 



Figure 2. P- hyperfine splitting results. 



charmonium 5-hyperfine splitting by at least 30- 
40%. Some have suggested that much of the 
remaining discrepancy is due to quenching er- 
rors. Our results indicate this is unlikely to be 
the case — we find at most a 10% effect from un- 
qucnching in the hyperfine splittings with both 
0(v 4 ) and 0(v 6 ) NRQCD. 

In contrast, the choice of tadpole factor and 
evolution equation each affect hyperfine splittings 
by 20 MeV or more — the instability of the evolu- 
tion equation should be investigated further. Fur- 
ther, 0(a s ) radiative corrections may shift the 
spin- splittings by as much as 40%. While this is 
a crude estimate, such sizeable corrections could 
easily swamp the small quenching error, and need 
to be carefully analysed. 

Thanks to Howard Trottier and Randy Lewis 
for stimulating discussions. 
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